Evidence of natural isotopic distribution from single-molecule SERS 
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We report on the observation of the natural isotopic spread of carbon from single-molecule Surface 
Enhanced Raman Spectroscopy (SM-SERS). By choosing a dye molecule with a very localized 
Raman active vibration in a cyano bond (C=N triple bond), we observe (in a SERS colloidal liquid) 
a small fraction of SM-SERS events where the frequency of the cyano mode is softened and in 
agreement with the effect of substituting ^^C by the next most abundant ^'^C isotope. This example 
adds another demonstration of single molecule sensitivity in SERS through isotopic editing which 
is done, in this case, not by artificial isotopic editing but rather by nature itself. It also highlights 
SERS as a unique spectroscopic tool, capable of detecting an isotopic change in one atom of a single 
molecule. 

PACS numbers: 



As a spectroscopic technique, single-molecule Surface 
Enhanced Raman Scattering (SM-SERS) is now well es- 
tablished and accepted. A precis of some of its most 
salient aspects has been recently provided in Ref. [T]. 
Arguably, one of the most refined and elegant demon- 
strations of Single Molecule Surface-Enhanced Raman 
Scattering (SM-SERS) is through the bi-analyte SERS 
technique [3 El Hj, using isotopically edited molecules 
[5l [6] . Isotopic editing is done purposely on specific moi- 
eties of standard SERS probes (like rhodamines [SI 'BJ ) to 
obtain molecules that have nominally identical chemical 
properties but distinguishable Raman features. 

Still, the typical atoms that constitute the structure 
of standard dyes used for SERS have their own natural 
isotopic spread. Accordingly, we can ask if the reverse 
logic applies: under SM-SERS conditions, is it possible 
to discern natural isotopologues? In this paper we show 
that -for specific cases- the natural isotopic spread of 
carbon in organic dyes is indeed detectable. 

One of the obvious contributions to the isotopic distri- 
bution of organic molecules (of interest for SERS) is car- 
bon. Despite a relatively small natural isotopic spread 
between (98.9%) and (1.1%) [7], carbon still 
plays a decisive role in the isotopic spread of organic 
molecules due to it preponderance with respect to other 
atomic species. Hydrogen is, of course, very abundant 
too in typical organic dyes, but this is compensated by its 
much smaller natural isotopic spread between hydrogen 
(99.985%) and deuterium (0.015%) 0. Let us study one 
particular example of a useful SERS probe: rhodamine 
800 (RII800). The reason for choosing this particular 
molecule will become clearer later. The structure of 
RII800 is shown in the inset of Fig. [l] it has the chemical 
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formula: C26H26N3O"'", and a molar mass of 396.21 g/mol 
when all the atoms are in their most abundant (lightest 
here) versions. Carbon and hydrogen are the most abun- 
dant atoms in the structure. Nitrogen and oxygen are not 
only a minority, but they also have smaller natural iso- 
topic spreads with respect to carbon (^^N: 99.63%; -^^N: 
0.37%, and ^^O: 99.75%; "O: 0.05%; ^^O: 0.20%) [7]. It 
is quite clear that by having the largest abundance and 
the widest (relative) natural isotopic spread (even when 
it is only '--^ 1%), it is enough for carbon to play the most 
predominant role. In fact, the "combinatorics" of iso- 
topic mass distributions are straightforward to calculate 
and can be compared to direct experimental measure- 
ments of the mass distribution of RH800 molecules, a 
common practice in mass spectroscopy analysis. This is 
explicitly shown in Fig.[l] obtained by standard high res- 
olution mass spectroscopy on our sample. Note that mass 
spectrometry typically uses large (macroscopic) amounts 
of molecules to discern these differences. 

These natural isotopic changes in mass result, in prin- 
ciple, in slight changes in vibrational frequencies. How- 
ever, these cannot be resolved in most cases for a variety 
of reasons that include: (z) the natural isotopic spread 
is very small, (ii) natural isotopic substitutions occur at 
random places within the molecular structure, and not 
necessarily at places where they will considerably affect 
the eigenvectors and frequencies of Raman active vibra- 
tions, and (iii) natural isotopic shifts are in many cases 
too small and hidden within the homogeneous linewidth 
of Raman peaks. In fact, a big fraction of the eigenvec- 
tors that produce Raman active modes in medium-size 
or large molecules will be extended over many atoms [S] , 
and the perturbation of a small mass change at a certain 
site will be relatively minor. From the standard theory 
of vibrations in molecules, the frequency of a given mode 
Ui in a molecule is related to its reduced mass fii through 
the proportionality uji cx (iJ-i)^^ [8 . The reduced mass 
can be viewed as a eigenvector-weighted mass-average 
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FIG. 1: Experimental (gray) and calculated (red) isotopo- 
logues for rhodamine 800 (RH800, inset) produced by the nat- 
ural isotopic spread of its constituent atoms. The experimen- 
tal results have been obtained by high-resolution mass spec- 
troscopy (MS). Approximately ~ 20% of RH800 molecules 
have one more unit of mass because of the isotopic substitu- 
tion of a ^^C^^^C in the structure (the rest corresponding to 
other isotopic substitutions). The number of molecules that 
have the ^^C^^^C substitution at the cyano bond (C=N) is, 
however, much smaller (~ 1%) because it is only one site; i.e. 
it only depends on the natural isotopic ratio of carbon. 




FIG. 2: (a) Molecular structure of RH800. In (6) we 
show the immediate environment of the cyano bond (C=N,) 
which forms a localized high frequency vibration resembling a 
"dumbbell" (inset). This vibration is Raman active, as shown 
in the spectrum in (c), and appears at ~ 2230 cm~^ in the so 
called "Raman silent" region (only triple bonds have Raman 
active vibrations in this region). Being a high frequency lo- 
calized vibration, the cyano stretching mode is particularly 
susceptible to the isotopic substitution ^^C^^'^C, resulting in 
a frequency softening of the mode by Au ~ 55cm~^. 



that takes into account the relative participation of dif- 
ferent atoms in the collective motion described by a par- 
ticular eigenvector. A change of mass Arrij on the j-th 
atom results in a relative change of frequency uji given 
by: Aujf — — Wi/(2/ii)A/i^ with A/i^ = {dfii/dmj)Amj. 
In general, this change will be small, and within the in- 
trinsic linewidth (homogeneous broadening) of the peak. 

There are nevertheless exceptions, a family of which 
is associated with the existence of localized vibrational 
modes (for which dfii/drrij is much larger). RH800 
has such a mode. The carbon- nitrogen triple bond 
(cyano group) in the structure of RH800, highlighted in 
Fig. [2]^ a), is a strong bond that produces a relatively iso- 
lated and localized stretch vibration. The cyano group 
(highlighted in Fig. ^h)) produces a Raman active vi- 
bration at a frequency of w ~ 2230 cm~^ for the most 
common isotopic combination of ^^C=^*N. This vibra- 
tion is fairly localized to the stretch motion of carbon 
against nitrogen and occurs in the so-called "Raman 
silent" region; above the standard fingerprint region (~ 
100 — 1700 cm~^) and below the hydrogen-stretching re- 
gion (~ 2900 — 3100 cm~^) for typical organic molecules. 
The cyano stretch Raman- active mode is shown explic- 
itly in the spectrum in Fig. [2j^c) for the specific case of 
RH800. Only triple bonds produce Raman active modes 
in this region. To a very good approximation, it is possi- 
ble to think of the cyano bond as a "dumbbell" formed by 
the C and N atoms, coupled (through a "softer" spring) 
to the rest of the molecular structure. This is shown 
explicitly in the inset of Fig. [2] By being a localized vi- 
bration with a high frequency, the C=N is particularly 



susceptible to isotopic substitution and, in particular, to 
the most common ^^C^^'^C replacement. The identifi- 
cation is helped by two facts: (i) the shift also occurs in 
a region where there cannot be any overlap with other 
Raman active modes, thus facilitating its identification, 
and (ii) the shift is large (because it is proportional to the 
frequency of the mode) and the peak is narrow enough to 
produce a new peak that cannot be confused as a small 
shift of the normal case (that could have other origins, 
like anharmonic/thermal effects). 
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TABLE I: Carbon-nitrogen (triple bond) stretch for different 
isotopic combinations in cyanobenzene calculated by Density 
Functional Theory (DFT). The most likely case observed ex- 
perimentally is the bond ^^C='^^N, followed by ~ 1% cases 
of ^^C=^^N which introduces a (calculated) frequency soften- 
ing of All) ~ 55cm~^. See the supplementary information for 
more details on the DFT calculations. 

Considering initially the most abundant case: 
i2(-<^i4]\j^ the reduced mass of the isolated "dumbbell" 
changes by ~ 0.28 for ^'^C^^'^G. . Hence, the 
predicted frequency change in the cyano frequency 
^c=N ~ 2230 cm~^ is Au;c=n ~ 50 cm^^. More accurate 
values that consider the link of the "dumbbell" to the rest 
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of the structure can be obtained by direct calculations 
with Density Functional Theory (DFT). One advantage 
of the the cyano bond is that, being a localized vibration, 
it can be modeled with any molecule that resembles the 
local chemical environment of the bond. For example, 
we can study the cyano bond dynamics in RH800 with 
calculations on the much smaller molecule cyanobenzene 
(CgHsCN) with obvious computational advantages. 
Further details of the calculations are provided in the 
supplementary information. A summary of the results is 
provided in Table |l] A frequency softening of ~ 55 cm^^ 
between the most abundant form of the bond (^^C=^^N) 
and next most abundant one (^'^C=^^N) is predicted by 
DFT. Such a shift is much larger than the linewidth of 
the corresponding Raman peak and therefore has to be 
easily observable in a small fraction of the SM-SERS 
spectra. 
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FIG. 3: SERS spectra of RH800 showing from top to bot- 
tom: (i) the normal spectroscopic feature of the Raman mode 
of the ^■^C=^'*N-cyano bond at ~ 2230 cm"'^ (top 2 spectra, 
red), (a) Mixed cases (next 3, black) where both the '^^C='^''N 
and ^^C=^*N-cyano bonds can be seen. These mixed spec- 
tra appear in the standard bi-analyte SERS method and we 
show three cases with larger, equivalent, or smaller relative 
intensities of one peak with respect to the other. (Hi) Single- 
molecule SERS events for the ^^C^^^N -cyano bond (next 2, 
blue). We also show the average spectrum (bottom, green) 
over the 18000 obtained spectra. Note that the number of 
cases where the ^^'C^^^'N -cyano bond signal is observed is 
very small, as evidenced by the negligible effect on the aver- 
age spectrum. We show the spectra in a range that includes 
the ~ 1650 cm~^ mode of RH800, which remains unchanged 
in all cases (as expected). 



SERS spectra have been collected in immersion for 
RH800 in a silver (citrate-reduced) Lee & Meisel colloid 



fS] at 10 mM KCl. The experimental conditions were 
identical to those reported in our recent study of iso- 
topologues of rhodamine in Ref. [B]. Spectra were ob- 
tained with a X 100 objective indexed-matched to water 
and the 633 nm line of a HeNe laser. Data were collected 
by a Jobin-Yvon LabRam spectrometer equipped with 
a notch filter and a nitrogen-cooled CCD. The dye con- 
centration in the colloidal solution was 10 nM and 18000 
spectra were taken with 0.2 sec integration time (with a 
1 sec dwell time in between spectra, to ensure statistical 
independence). We know from previous characterizations 
of our system for similar experiments with the bi-analyte 
method that we are in the conditions where many of 
the spectra are single molecule in nature, with a frac- 
tion of them having contributions from more than one 
molecule. 

Once the spectra are taken, we need to find cases where 
the isotopically modified cyano bond is present. This 
implies looking for a "minority" of cases that will not 
contribute much to the average signal. We have done 
this using Principal Component Analysis (PCA) (as ex- 
plained in the supplementary information), but this is 
not necessary and any method should render equivalent 
results. The big advantage of PCA is that it reduces the 
search for the appropriate spectra to a matter of seconds; 
a task that becomes very tedious otherwise when 18000 
spectra need to be analyzed. 

Figure [3] shows representative spectra of the different 
cases of interest here. The vast majority of events show 
the cyano bond Raman peak at ^ 2230 cm~^. A closer 
look, however, reveals situations with mixed signals at 
the corresponding frequencies expected for ^^C=^'*N and 
^^C=^'*N. These are equivalent to the mixed signals that 
appear in the standard bi-analyte method. The exper- 
imental frequency softening is ~ 56cm~^, in excellent 
agreement with the DFT predictions (Table |l]). We chose 
three spectra in Fig. [3] showing a larger, equivalent, or 
smaller intensity of the the -'^^C=^^N-peak with respect 
to that for ^'^C="'^*N. Finally, we show cases where the 
signal comes only from the isotopically modified cyano 
bond (^-^C^^^N). The average SERS spectrum is also 
shown for direct comparison. It is interesting to note 
the small feature at ~ 2120 cm~^ with an integrated in- 
tensity « 10 times less than the cyano bond peak. This is 
assigned to an overtone band and has no relation to the 
isotopic spread. It does however give an idea of how small 
the ^■^C=^'*N Raman peak should be in the average sig- 



nal: 



10 times smaller than this ~ 2120 cm overtone 



band, i.e. clearly undetectable in the average spectrum 
under the present conditions. This is in stark contrast 
with the SM-SERS spectra, where the (i3CEE"N)-cyano 
bond signal can completely dominate the spectra. The 
single-molecule nature of the signals here enable us to ac- 
cess information, such as the isotope-induced shift, that 
are clearly washed out in the ensemble average. 

The results of this paper provide inderectly another 
demonstration of SM-SERS sensitivity to the ones al- 
ready known [Tl [TUl |TT]. But with this aside, the main 
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claim here is one step beyond another demonstration of 
SM-SERS sensitivity. We have indeed shown a specific 
example of how single-molecule SERS can be used as tool 
for the observation of weak spectroscopic features that 
are otherwise washed out in ensemble averaged spectra. 
Other types of spectroscopy (like fluorescence) would be 
completely insensitive to a change of one unit in the iso- 
topic mass of one atom, even if single molecule sensitivity 
is achieved. We believe the observation of naturally oc- 
curring single-molecule isotopic fluctuations in SERS is 
a pleasing demonstration of the consistency of SM-SERS 
phenomena in revealing subtle aspects of molecular spec- 
troscopy that cannot be observed otherwise. The spectra 
in Fig. [3] show, effectively, examples of single-molecule 
Raman spectra of molecules that differ from the others 



by one unit of atomic mass in one atom. Finally, it is also 
worth pointing out that the observation of natural iso- 
topic spread with SM-SERS is not necessarily restricted 
to triple bonds. Options with chlorine-containing dyes 
could be particularly interesting to explore, for the natu- 
ral isotopic spread of chlorine can produce a much more 
balanced list of naturally occurring isotopologues. 

We are indebted to Sean Buchanan for help with the 
measurements. PGE and ECLR acknowledge partial 
support from the Royal Society of New Zealand (RSNZ) 
through a Marsden Grant. 

Supplementary information available: details of DFT 
calculation and PCA analysis of the SERS spectra. 
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